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Abstract 

We show how to couple two critical Q-state Potts models to yield a new self-dual critical 
point. We also present strong evidence of a dense critical phase near this critical point when 
the Potts models are defined in their completely packed loop representations. In the continuum 
limit, the new critical point is described by an SU(2) coset conformal field theory, while in 
this limit of the the critical phase, the two loop models decouple. Using a combination of 
exact results and numerics, we also obtain the phase diagram in the presence of vacancies. 
We generalize these results to coupling two Potts models at different Q. 

1 Introduction 

The Q-state Potts model is of fundamental interest in two-dimensional statistical physics. It 
arises in an amazing number of different contexts in mathematics and theoretical physics alike: 
algebra, combinatorics, probability theory, graph theory, integrability, conformal field theory 
(CFT), stochastic Loewner evolution (SLE), etc. Thus progress in understanding the Potts 
model serves as a useful gauge for progress in many of these fields. 

Via its underlying algebraic formulation — in terms of the Temperley-Lieb algebra [1] — the 
Potts model can be recast in several similar ways, but which do however present subtle differences. 
Its formulation as a model of completely packed loops with non-local weights is prototypical within 
the large (and still growing) class of so-called loop models in classical statistical mechanics. Such 
models push the limits of CFT, as they present non- unitary and logarithmic features [2j, and 
they link directly to physics-inspired models of extended objects, such as polymers. The loops 
underlying the Potts model role as convenient lattice regularizations of SLE traces [3] . They can 
also be interpreted as particle trajectories in various models of disordered systems, such as the 
spin quantum Hall effect [4]. 

Studying several coupled Potts models is a challenge for all of the techniques listed above, and 
also serves as a worthwhile purpose on its own right. 

In particular, coupled Potts models arise naturally when studying quantum loop models [5j[6]. 
Each configuration in the classical loop model is a basis element of the Hilbert space of the 
quantum loop model. The reason why two coupled copies of the loop model arise in this context 
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is elementary quantum mechanics. When computing expectation values, one of course weights by 
the wavefunction squared. Thus one loop model comes from the bra, and the other from the ket. 
Our results are very useful for quantum loop models, because the critical phase we find below 
in the classical model implies the existence of a phase with topological order in the quantum 
loop model with completely packed loops [B]- In a topological phase, the excitations are 
fractionalized, and have anyonic statistics. These features are possible when the ground state of 
the quantum system is a sum over different loop configurations; the non-locality of loops makes 
anyonic braiding possible. For the anyonic excitations to be deconfined, the underlying classical 
loop model must be critical, so that all lengths of loops appear in the partition function. 

Coupled Potts models are also at play when studying a single Potts model with quenched 
random bond disorder. Using the replica trick, this corresponds to coupling J\f pure Potts models 
via the product of their energy operators, followed by the formal limit N — > 0. For Q — 2 <C 
1 and in the limit of weak non-frustrated disorder, this can be done within the perturbative 
renormalization group (RG) approach of CFT |8j. There exists ample numerical evidence that 
the replica symmetry is not broken in this case. Moreover, the formal expressions obtained for 
the critical exponents corresponding to N > 2 integer have been identified — order by order in 
the perturbation — with those of particular critical points of N coupled lattice models. If exact 
results for these critical points could be obtained exactly, their analytical continuation to the 
M — > limit would provide exact results for the random-bond Potts model [9]. In this context, 
the case of two coupled models to be studied here is somewhat particular, since the RG flows 
depend on the coupling g in the combination (2 — N)g. 

In the present study we assume (as in [9]) that the coupled Potts models are self-dual, i.e., 
they are left invariant under the combined duality transformation of both underlying models. 
Although it is known in the special case of Q = 2 that mutually dual pairs of critical points do 
exist, it seems plausible that self-duality will generically allow for the highest multicritical — and 
thus most interesting — behavior. 

Among the directions that we leave for future investigation, one of the most elusive is un- 
doubtedly to provide a Coulomb Gas (CG) description for coupled Potts models. A putative CG 
construction would allow one to compute the critical exponents of certain non-local operators 
analytically, an issue that we only pursue numerically here. Suffice it to say that our analyt- 
ical result for the central charge at the new self-dual critical point hints at a CG with three 
height components, but only two of those are easy to identify geometrically. This is not without 
reminding us of Ref . [10] . 

The plan of the paper is as follows. In section [2] we introduce the models to be studied, and 
summarize our main results. In section [3] we derive the location and universality class of a new 
self-dual critical point. The cornerstone of the arguments is a level-rank duality in a related 
SO(N) symmetric integrable model. Numerical calculations are used to verify the consistency of 
the proposed phase diagram, and to give results for the dimensions of watermelon operators at the 
new critical point. In section [J] we argue that adjacent to the new critical point, the loop model 
possesses a critical phase, with parameter-independent critical exponents. We provide evidence 
that the phase transition into this phase is simultaneously first and second order. In section [5] we 
generalize the model to include vacancies. We argue that dilution is an irrelevant perturbation of 
the completely packed model. The phase diagram for the dilute model is presented. Finally, in 
section [6] we generalize our results to the case where the two coupled Potts model have a different 
number of states. 
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2 The models and their phase diagrams 



The original formulation of the Q-state Potts model places a spin <Xj taking integer values 
1,2, ... ,Q on each site of some lattice. In this paper we shall discuss exclusively the square 
lattice, but we expect by universality that our results will apply qualitatively and sometimes 
quantitatively to other lattices. The Potts model has nearest-neighbor interactions invariant un- 
der the permutation group Sq, so that the interaction energy J8 ait(T . only depends on whether 
adjacent spins (Xj and a a are the same or different. 

The Potts model on the square lattice has a Kramers- Wannier (high- to low-temperature) 
duality This is a non-local map of spins which preserves the bulk free energy. The self-dual 
point is critical [12] when Q < 4. In a field-theory description of the region around these critical 
points, taking the temperature off the self-dual value corresponds to perturbing by the energy 
operator e (the continuum limit of J5 (7u(7j ). The energy operator is odd under duality, e — > —e, 
so that perturbing by this operator with one sign describes the high-temperature phase, while 
the other sign gives the low-temperature phase. 

Although the original definition of the Potts model requires Q to be an integer, there exist 
generalizations to all Q which are still commonly known as the Potts model. We will treat two 
closely-related but distinct models, which we refer to as the loop and the height models. The 
transfer matrices of the original Potts model and the two generalizations are most conveniently 
defined in terms of the Temperley-Lieb algebra; each model then corresponds to a different 
representation. The algebra for a system of width L has L generators e« acting at positions 
i = 1, 2, . . . , L as well as the identity /, which obey the relations pQ 

ef = \fQei, e%ei±\ei = e iy eiej = ejCi for \i - j\ > 1 (1) 

The transfer matrix at the self-dual isotropic point for width L is then 

T = R1R3 . . . RL-1R2R4 ■ ■ ■ Rl (2) 

where Ri = I + et for the ferromagnetic self-dual point, and R = I — for the (non-unitary) 
antiferromagnetic self-dual point. The partition function for an M x L lattice with periodic 
boundary conditions in the "time" direction is then given by 

Z = tr T M . (3) 

In this formulation, duality amounts to exchanging with / on each plaquette. 
It is convenient to parameterize 

defining the parameter k. When k is an integer obeying k > 2 (so that 2 < Q < 4) representations 
of the algebra have special properties. One important one discussed below is the existence of a 
height representation, with positive and locally defined Boltzmann weights. Having k integer is 
also desirable for the quantum loop models, because the existence of the Jones- Wenzl projector 
at these values means that the number of ground states remains finite even when the theory is 
defined on a torus [7] . 

In the loop representation of the Temperley-Lieb algebra, the generators act on L strands, 
as shown in Fig. [TJ The result is a completely packed set of self-avoiding and mutually-avoiding 



3 



i9 <8> 



Figure 1: Action of I (left picture) and e% (right picture) on a pair of strands at positions i and 
i + 1. The transfer (time) direction is upwards. Each picture corresponds to a plaquette of the 
dual lattice. 

loops. Complete packing means that the loops cover every link of the square lattice; the only 
degrees of freedom are the two ways in which the loops avoid each other at each vertex of the 
lattice. The transfer matrix adds L plaquettes to the dual lattice as illustrated in Fig. [21 so it 
is acting on a zig-zag row of L plaquettes. Multiplying out the product in ([2]) makes it obvious 
that each configuration in the completely packed loop model appears with a coefficient which is 
just (VQ) #loops when R4 = I + e», and (-^/Q) #loops when B4 = I - e*. Note in particular that 
the parameter k in Q does not need to be an integer for the loop representation to be defined. 




Figure 2: Action of a term in T 2 on a zig-zag row of L plaquettes. The left (resp. right) picture 
illustrates periodic (resp. free, or reflecting) boundary conditions in the space direction. 

A remarkable property of the loop representation is that the Boltzmann weight of a config- 
uration depends only on topological properties of the loops. (In fact, for this reason the loop 
representation of the Temper ley-Lieb algebra underlies the Jones polynomial in knot theory.) 
The local topological properties follow directly from the Temperley-Lieb algebra ([I]). The second 
relation means that configurations which are topologically equivalent receive the same weight, 
while he first relation means that a closed loop not wrapping around a cycle receives a weight 
y/Q relative to the configuration without the loop. 

To evaluate the full partition function in the loop representation, we need to define the trace 
in ([3]), which corresponds to imposing periodic boundary conditions in the "time" direction. It 
is natural to require that closed loops wrapping around the time direction also receive a weight 
\[Q. Thus the weight of a given loop configuration C is found simply by counting the number 
ri£ of loops is the configuration: 

z = Y J Q ncl2 
c 

This defines the Markov trace of elements of the Temperley-Lieb algebra. A loop is inherently 
a non-local object, so even though the Boltzmann weights here are all positive for Q > 0, the 
corresponding field theory is typically not unitary. 

When k is an integer, the Temperley-Lieb algebra has a height representation where the 
resulting lattice model (often known as RSOS, or restricted solid-on-solid, model) has positive 
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Figure 3: The four configurations at each vertex in the coupled completely packed loop model. 



Boltzmann weights |13|, [T^ [T5l 116] . The heights h{ = 1, 2, . . . , k + 1 live on the dual lattice, so 
that each plaquette illustrated in Fig. [5] has four heights at its corners. Heights of neighboring 
sites satisfy the RSOS constraint — hj\ = 1. The generator ej now acts locally as 
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,hi,hi+i) = S hi _ uhi+1 2^ ■ / ; n I \hi-x,K>hi+i) ( 5 ) 




and the trace in ([3]) is simply the usual matrix trace. 

The connection between the loop and height models — valid when both are defined, i.e., for k 
integer — is based on exact partition function identities [15J that hold strictly speaking only when 
the underlying lattice is planar. Imposing doubly periodic boundary conditions leads to further 
subtleties, which are however not important in the context of this work. 

This paper is devoted to studying two coupled self-dual Potts models, in both their loop 
and height formulations. We thus consider two Potts models on the same lattice, described 
respectively by two independent sets of Temperley-Lieb generators labeled by e, and /j. Here 
we require that Q be the same for the two models; later we will consider the more general 
case. We require that the coupled model remain self-dual under the combination of the duality 
transformations in both models, and that it be invariant under exchange of the two models. The 
transfer matrix remains of the form ([2]) , with Ri for each plaquette given by 

Ri = I <g> I + <a <g> fi + A (I <8 fi + e; <g> I) (6) 

In the loop representation, this is very simple to represent; the four possible configurations 
at each vertex are displayed in Fig. [3l The two copies mean that on each link there are two 
species of loops. Loops of each species cannot cross at vertices, but the two species can cross. 
The partition function is a sum over all configurations C and M. of the two species of loops: 

Z = J2Y1 Q {nc+nM)/2 X nx (7) 
C M 

where nx is the number of times the rightmost two vertices in Fig. [3] occur. This is a completely 
packed version of the loop model discussed in [T7] ; we will make contact with the results discussed 
there when we discuss the dilute generalization in section [5j 

When the parameter A = 1 and A = — 1 , we have decoupled ferromagnetic and antiferro- 
magnetic models respectively, both critical. When A = 0, Ri describes a single Q 2 -state Potts 
model, because the products ej ® fi generate a Temperley-Lieb algebra ([IJ with Q replaced by 
Q 2 . In the loop representation, A = amounts to requiring that the loops in the two models be 
identical, so that each loop gets a weig ht (VQ) 2 = Q. Taking the dual of one of the two models 
exchanges A i— > 1/ A, so we restrict our attention to |A| < 1. 
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In the height representation, the theory remains unitary for A > —\j\[Q. This follows by 
explicitly examining the Boltzmann weights; they all remain positive except for one type of 
plaquette. Namely, a negative Boltzmann weight occurs for A < for all plaquettes where the 
heights in one RSOS model differ in the time direction, while the heights in the other RSOS model 
differ in the space direction. Because of the periodic boundary conditions in the time direction, 
there must be an even number of such plaquettes. Thus one can flip the sign of these weights 
without changing Z, giving a local model with all positive Boltzmann weights. This remains true 
until A = —1/y/Q, where other weights change sign. 

For k = 2, perturbing A away from 1 is marginal and preserves the criticality. This is easy 
to see in the height formulation, because k = 2 corresponds to Q = 2: ([6]) then describes two 
coupled Ising models. This is known as the Ashkin- Teller model along its self-dual line, and the 
physics is well understood in the height formulation [18]. The model remains critical until A = 0, 
where it is equivalent to the critical 4-state Potts model. For A < 0, it is no longer critical. 

At the decoupled points for Q < 4, the models are critical and so described by a conformal 
field theory. In the height model for A = 1, the conformal field theories describing these critical 
models are known as the unitary minimal models, and have central charge 



'-^- jk+W+T) ) forA = 1 (8) 

The behavior of the decoupled height models at the antiferromagnetic critical point A = — 1 is 
somewhat trickier to describe. The negative Boltzmann weights make the theory non-unitary. 
However, if one instead uses a Hamiltonian formulation of a 1+1 dimensional quantum theory, 
the Hamiltonian corresponding to A = —1 is simply H = Yli=i e ii which is hermitian. The 
unitary conformal field theory describing the continuum limit of this Hamiltonian is known as 
the Zk parafermion model, whose central charge is 

Chei g ht = 2( 2( fc fc ~ 2 1) ) forA = -l. (9) 

In the loop model, (JSj) remains true for A = 1. However, since k is now a continuous parameter, 
the cases A = 1 and A = — 1 are now linked by analytic continuation, k i— > — ^qpp Note that under 
this transformation \[Q \— ► —yQ, which is tantamount to Ri = I + ej i— ► I — et as required. In 
particular, for A = — 1 the correct central charge is then obtained by analytically continuing ([8]) 
and reads 

6(fc + l) 2 ~ 



cioop = 2 ^1 - y k + 2 J \ for A = -1. (10) 

When k is integer, the corresponding critical theory completely disappears from Z due to massive 
cancellations [19] of eigenvalues in the Markov trace ([3]), and what remains is tantamount to the 
height model. 

Taking A away from 1 is a perturbation odd under each the individual dualities of the Potts 
models, and even under the combination. The simplest operator with these properties is given 
by the product of the energy operators £1^2 of the two models. This operator has dimension 
(k + 4)/(A; + 1), so it is relevant for k > 2; it is the only relevant one with the appropriate 
symmetries. Thus the question becomes whether this perturbation causes a flow to a non-trivial 
critical point. Compelling evidence suggests that there is no such critical point, although no 
definitive argument rules it out, Namely, this perturbation is integrable in the field theory [21] . 
and the only consistent S matrices describing this perturbed theory are massive, indicating there 
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Figure 4: Phase diagram for the loop model (top) and the height model (bottom). The arrows 
indicate renormalization group flows. The cross-hatching between the new critical point C and 
the decoupled antiferromagnetic critical point Da represents the critical phase. 

is no flow to a critical point. Moreover, this is consistent with the fact that at A = 0, the model 
for k > 2 is non-critical, because the Q 2 -state Potts model at its self-dual point is critical only 
for Q < 2. Thus it is natural to guess that for k > 2, the flow goes from A = 1 to A = 0, the 
latter being a stable non-critical fixed point of the renormalization group. Our numerical work 
supports this picture. 

The interesting new physics occurs for A negative, where we find a new critical point. We 
present this calculation in section [31 where we derive the exact location of the critical point and 
determine which conformal field theory describes its continuum limit. 

The arguments given above — in conjunction with numerical calculations presented in sec- 
tion 13.21 below — lead to the phase diagrams shown in Fig. SJ The point S is the strongly coupled 
point A = 0, which is non-critical for k > 2. The point C is the new critical point at A c , to be 
discussed in detail in section [3] below. The points Dp and Da are the decoupled ferromagnetic 
and antiferromagnetic critical points at A = 1 and A = — 1. The characteristics of Da differ 
between the height and the loop models. Finally, the point N is a non-critical point which occurs 
only in the height model, and which is required for the consistency of its RG flow diagram. 

Note in particular that the phase AG [—1, A c ) is non-critical in the height model, but critical 
in the loop model. The latter is an example of a "Berker-Kadanoff " phase |19j . 

The point A = —l/y/Q, the lowest value of A where the height representation is unitary, is 
special. It is quite important for quantum loop models, and it lies in the critical phase for k < 6, 
and is the new critical point for k = 6. When A = the weight of each configuration 

can explicitly be given in terms of a chromatic polynomial squared [6\. This can be seen for the 
square lattice by noting that 



where d = y/Q- The projectors 1 — e/d and e/d are orthogonal. We divide the square lattice 
on which the loops live into the usual two sublattices (even and odd). The Zou>-temperature 
expansion of a single Potts model is found from writing i?p tts = (1 — e/d) + (d — l)e/d on 
sublattice 1, and i?p tts = (e — 1/d) + (d — l)/d on sublattice 2, and then associating a domain 
wall with each 1 — e/d term included on sublattice 1, and e — 1/d on sublattice 2. It is easy 
to check that these walls cannot end, although there can be three or four walls touching each 

x We note, however, that the S matrices given in [5T] for k > 2 are wrong — they are the level-rank duals of the 
correct ones, in the technical sense described below. 



R = (1 

= (c 



e/d)(l - f/d) + (d 2 - l)(e/d)(f/d) 
l/d)(f-l/d) + (d 2 -l)(l/d 2 ) 
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site. If we sum over the spins for any fixed domain-wall configuration, the weight per graph is 
proportional to the chromatic polynomial of the dual lattice. This is easiest to see by simply 
proving it for all Q integer in the original Potts representation; by definition of the chromatic 
polynomial it must apply for all Q. Thus coupling Potts models at A = —1/d corresponds to 
requiring that the low-temperature expansion of the two be the same. (This is proved for any 
lattice with four links touching each site in [6].) The behavior of the low-temperature expansion 
at A = —1/d is therefore analogous to that of the completely packed (Fortuin-Kasteleyn [2D]) 
loop expansion at A = 0. 



3 The new self-dual critical point 

In this section we show that at a particular value A = A c , the coupled Potts models have a critical 
point. In accord with the discussion of the previous section, this value must be negative, and we 
show here that 

1111 



(12) 



The conformal field theory describing this critical theory is the coset model [22 

SU(2) k x SU(2) k 
SU(2) 2k 

with central charge 



3A- 
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(fc + !)(*; + 2) 
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3.1 Level-rank duality and the BMW algebra 

At the value A = A c , the coupled theory is integrable. We show this by using level-rank duality to 
write the transfer matrix here in terms of another algebra, the SO(N) Birman- Murakami- Wenzl 
(BMW) algebra [23]. Lattice models involving this algebra were introduced long ago [23], and it 
is possible to show that integrability follows directly from algebraic properties, and so applies to 
any representation [25]. In this more general setting, the Andrews-Baxter- Forrester height models 
[T3] discussed are associated with SU{2) and its quantum-group algebra deformation U q [SU(2) s j . 

The BMW algebra was introduced to study representations of the braid group and find new 
knot and link invariants. Braid generators in general satisfy BiBj = BjBi for \i — j\ > 1, and 

BiB i+ iBi = B i+ iBiB i+ i . (14) 

This is known as the "third Reidemeister move" in knot theory. The braid generators can be 
inverted: B{B^ = / (the second Reidemeister move). A useful graphical representation of 
the braid group is given by each Bi and B^ 1 as acting on two strands like the Temper ley-Lieb 
generators, but where B and B^ 1 correspond to overcrossings and undercrossings respectively. 

The SO(N) BMW algebra is generated by the B { and Br 1 , where i = 1,2,..., L. This 
algebra involves a parameter q as well as N; for the models we label SO(N) k , q = e t7T /( N + k ~ 2 ) 
is a root of unity. The braid generators of the BMW algebra of course satisfy the second and 
third Reidemeister moves, but the algebra has more structure. The Temperley-Lieb algebra is a 
subalgebra, with Ei defined as 

Bi -Br 1 



where q is a parameter. Here the Ei obey 



{Eif = DEi, EiE i+1 Ei = E h E { Ej = EjEi for \i - j\ > 1 

where 

9-9 

The full set of BMW relations can be found in |25|. One useful one is the first Reidemeister move 



BiEi = EiBi = q N x Ei . 

For later use, it is convenient to define 

Xi = q~ l I + qEi - Bi (16) 

This generator Xi is particularly interesting in the N = 3 case, where the BMW algebra is related 
to the chromatic algebra [26j EZ] • 

The coupled Potts model transfer matrix © can be written in terms of the £0(4)^ BMW 
algebra. This BMW algebra is special because the algebra SO (A) is the direct product of two 
SU(2) algebras. One should think of the strands in SO(N)k as behaving in the vector repre- 
sentation of the quantum-group algebra U q (SO(N)). As with ordinary Lie algebras, U q (SO(A)) 
decomposes into copies of U q (SU(2)), so the lines in 50(4)^ behave as two independent spin- 
1/2 lines in SU{2)^. Spin-1/2 lines obeying the Temperley-Lieb algebra satisfy the braid group 
relation 6j6j + i6j = bi + ibibi + i and bib^ 1 = I if we define the bi via the "skein relation" 

h = <T 1/2 - 9 1/2 e*, &7 1 = 9 V2 - Q- 1/2 e h 

where the ej obey ([I]) for yfQ = q + g -1 . (This is the braid-group representation underlying the 
Jones polynomial.) So less abstractly, given two independent representations labeled and fi of 
the TL algebra, we can construct one of the SO (4) & BMW algebra by defining 

E^ = ei ®fi (17) 
B^ = (q-Wl - q^e t ) (q-^I - q^ /<) = q^I + qE® - X® (18) 

where given the definition of Xi in (|16fl . we have 

Xf =ei®I + I®fi 

for this representation of #0(4)^. Therefore the transfer matrix ([6|) for each plaquette can be 
written simply as 

R l = I + Ef + \X\ A) . (19) 

Because of the relation (I15p . we need use only three of the generators for each plaquette. 

Lattice models whose transfer matrices are built on the SO{N)k BMW algebra are integrable 
only for special values of the parameters. Height models associated with all non-exceptional Lie 
algebras were introduced in [23]. The transfer matrices for the critical points of these lattice 
models can be written in terms of BMW algebras; for the SO(N)k case they are [25J 

R<i N) («) = [7 - u] [1 - u]I + [u] [1 - 7 + u]Ei + [ 7 - u] [u]Xi (20) 
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where 7 = (JV — 2) /2 and we have introduced the shorthand 



The parameter u is simply a lattice anisotropy which does not affect any of the universal prop- 
erties. Under 90 degree rotations, I and E are interchanged, while X remains invariant. Thus 
isotropic Boltzmann weights are found by setting u = 7/2, giving 

*r'(7/2) _ I+E . m Xi (21) 



[ 7 /2] [l- 7 /2] [I-7/2] 

Prom the BMW algebra it follows that this R matrix satisfies the Yang-Baxter equation 

Ri(u)R i+ i(u + u')Ri(u) = R i+ i(u)Ri(u + u')R i+ i(u) , 

so that the model is integrable. Note that the braid generators are recovered in the u — ► ±00 
limits: 

Ft {N) (u) 

hm %^= q ^I + q ^E l -X l = Bf . 
u— >±oo q^\ u 1) 

It is straightforward to establish that lattice models with the transfer matrix given by (|20p 
and (|2|) are critical. This is a generic feature of integrable models with trigonometric Boltzmann 
weights. Here one can see this by utilizing a deformation of these models which preserves the 
integrability [2jj . This enables one for example to compute the free energy by using the inversion- 
relation method [18] . From this one sees that as the deformation parameter is tuned back to the 
trigonometric weights, the free energy behaves non-analytically, as it does at a critical point. The 
central charge of the conformal field theory describing this critical point can be computed after 
making some (very believable) technical assumptions; for u > it is [28] 

N ( N 2 — 3N + 2 

c = — 1 



2 V (N + k-2)(N + k-3) 
which is the central charge of the 

SO(N) k ^ x 5Q(JV)! 

SO(N) k { > 

coset models. 

The new critical point for coupled Potts models does not correspond to the N = 4 case of 
(|21|) . Rather, this yields the decoupled ferromagnetic critical point. Namely, when N = 4, 7 = 1, 
and the decomposition of the SO (4) BMW algebra in terms of two Temper ley-Lieb algebras in 
(I17I18II means 

Rf\u) = ([1 - u]I + [u] ei ) ® ([1 - u]I + [u]fi) 

Thus the models decouple here: at the isotropic point u > this is simply the A = 1 ferromagnetic 
critical point. Indeed, because as conformal field theories 50(4)^ = SU{2) k x SU{2)j t} the coset 
in (I22p simply splits into two copies of 



SU(2) k -i x SU(2)i 
SU(2) k 
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with central charge ([8]). 

To find the new critical point at A = A c , we need to exploit the level-rank duality of the BMW 
algebra [27]. It is simple to check that if the Bi satisfy this SO(N) k algebra, then the —B^ 1 
satisfy the SO(k)j^ algebra. This follows because q is invariant under this exchange, and 

so that D remains unchanged. Thus from any representation of SO(N) k with generators {Bi, Ei}, 
we can construct one of SO(k)N with generators {Bi,Ei} by 

Bi = —B i , Ei = Ei. 

We then have 

g- 1 / + qEi + Br 1 

q-H + qEi + {q- q^I - E t ) + B, t 
{q + q' l ){I + Ei)-Xi 

The transfer matrix of the coupled Potts models defined in ([6|) has been expressed in terms 
of the SO{A) k BMW algebra via (|19j) . The level-rank duality therefore implies that it can also 
be written in terms of the generators {Ei,X{\ of the SO(k)4 BMW algebra: 

Ri = {l + \{q + q- 1 )){I + Ei)-\X i (23) 

We can now find the new critical point, because the SO{k)4 : models are critical when the Boltz- 
mann weights obey (|2ip . It follows that the coupled Potts models are critical at A = A c , where 
Ri oc Rf\ This holds if 

Ac [7/2] 
l + X c (q + q- 1 ) [1-7/2]' 

Plug ging in the appropriate values of [7/2] and [1 — 7/2] for SO{k) 4 and solving for A c gives (llip . 

We expect the central charge to be independent of the representation of the BMW algebra 
used, as long as appropriate boundary conditions are chosen. This critical point at A = A c 
therefore should have the central charge of the coset conformal field theory 

SO{k) 3 x 5Q(fc)i 
SO{k) 4 ' 

which is that given in (|13p . A related type of level-rank duality arises in coset conformal field 
theories [29] . This lets us map this coset model on to others: 

SO(k) 3 x 5Q(fc)i ^ SO(4) k _ SU{2) k x SU(2) k 
SO(k) 4 ~ 50(3) fc SU(2) 2k 

The latter coset is the one we gave above in (|12|) : note that the coset description involving SO (4) 
is not the SO (4) case of ()22|) . Thus exploiting the integrability of this lattice model lets us not 
only determine the location of the critical point exactly, but it lets us determine the correct coset 
conformal field theory. 

Note that when N = k, the level-rank duality does not change the model, only the couplings. 
The points A = 1 and A = A c are not self-dual under level-rank duality, but rather are mapped to 
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Figure 5: Central charge of the loop model with k = 3, as a function of A. The pink cross 
indicates the expected infinite-size behavior at the new self-dual critical point. 

each other. Thus for the coupled Potts models with N = k = 4 (Q = 3), the new critical point 
at A = A c is equivalent to the usual A = 1 pair of decoupled models. Indeed, the central charge 
(|13p here at A = A c is 8/5, that of two decoupled 3-state Potts models at their ferromagnetic 
critical points. As pointed out in section 5.2.2 of [9], the coupled Q = 3 models in fact have an 
exact symmetry J i— > — J, where J is the original Potts spin coupling. In the present notation 
we have A = (e J — l)/y / Q, and the symmetry reads A i— > — 1+ ^^q ■ For Q = 3 this gives indeed 

3.2 Numerical calculations 

The coupled model can be studied numerically using techniques for exactly diagonalizing the 
transfer matrix T, either in the loop or the height representation. We focus here on the loop 
model. 

For reasons of numerical efficiency, it turns out to be most convenient to chose the transfer 
direction to be axial with respect to the square lattice on which the loops live (as opposed to 
the diagonal transfer direction depicted in Figs. [2H3]). Periodic boundary conditions are imposed 
in the space direction. To avoid parity effects in the antiferromagnetic region, the number of 
strands L is chosen to be even. We can then handle widths up to L = 18. 

The central charge can be extracted from finite-size scaling of the leading eigenvalue of T 
in the standard way. The resulting three-point fits c(L — 4, L — 2, L) are shown in Fig. [5] as 
functions of A. The figure shows the case k = 3, but larger integer values of k lead to similar 
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Exponent 


k = 3 


k = 4 


k = 5 


£2,0 


0.998(1) 


0.999(1) 


0.999(1) 


2:2,2 


0.200(2) 


0.130(2) 


0.095(2) 


£4,0 


1.560(2) 


2.5 (2) 


— 


£4,2 


0.774(1) 


0.510(5) 


0.392(4) 


£4,4 


0.926(4) 


0.677(2) 


0.555(5) 


£6,0 


2.40 (3) 


2.47 (1) 


3.0 (1) 


£6,2 


1.70 (5) 


1.15 (5) 


0.84 (5) 


£6,4 


1.85 (5) 


1.28 (5) 


1.00 (5) 


£6,6 


2.08 (7) 


1.54 (5) 


1.26 (6) 



Table 1: Watermelon exponents x^i 2 at the critical point A c for k = 3,4,5. The error bar on 
the last quoted digit is given in parentheses. One entry is missing due to convergence problems 
in the diagonalization method being used. 

results. Eq. (jHJ) then predicts c = | = 1.4 at the decoupled ferromagnetic point A = 1, and 
Eq. (HU) gives c — — 1.35 at the new self-dual critical point, which according to (lllj) is 
situated at A c = — \/2sin(^) ~ —0.2212. Finally, for the strongly coupled point A = we have a 

single Q 2 = (2cos(^)) ~ 6.8541 state Potts model, which is non-critical but with a rather large 
correlation length. 

All of this is nicely brought out by Fig. [5j The fact that the peaks in c at the critical points 
A = A c and A = 1 narrow upon increasing L, further indicates that both points are RG repulsive. 

To summarize this far, our numerical results not only confirm the existence of the critical 
point at A = A c , but also provide compelling evidence that the central charge (I13p is the correct 
one. 

We now turn to the operator content of the critical theory. One of the most interesting 
and geometrically appealing operators that can be defined in a loop model is the so-called l-leg 
watermelon operator On [30J. To define this operator properly, we first imagine assigning an 
arbitrary orientation (clockwise or counterclockwise) to each loop in the model; the original loop 
model is then recovered by independently summing over the two orientations of each loop. The 
oriented loop model satisfies a zero-divergence constraint: the net arrow flow out of any simply 
connected region vanishes. The operator On breaks this constraint, and corresponds to a point 
(or a small neighborhood of a point) that acts as a source of I oriented loop segments. The critical 
exponents X£ corresponding to the asymptotic decay of the two-point correlation functions 

(Oi{xi)0- t (x 2 )) ~ |xi - x 2 |- 2 ^ for |xi - x 2 | « 1 (24) 

can be computed analytically when the loop model admits a Coulomb gas construction. This has 
been accomplished for the loop model corresponding to a single Potts model [30J as well as for 
more complicated loop models [10] . 

The context of J\f coupled loop models allows for generalizing these operators to {Ix^i^t ■ ■ ■ > ^v)~ 
leg watermelon operators ^ .^(x), where the pth. model supports an £ p -leg defect at the 
point x. The interpretation of this generalization in the replica setup was given in [9j. The 
corresponding critical exponents can be measured numerically by defining a modified transfer 
matrix in which the pth model is constrained to having precisely \£ p \ loops which wind around 
the time direction. 
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Returning to the case M = 2, we have determined numerically the exponents X£ lt e 2 (with 
£ p = 0,2,4,6 for p = 1,2) at the critical point (lllh . This study involved strips of width up to 
L = 16 strands. Results for parameter values k = 3, 4, 5 are shown in Table [TJ These results give 
very strong support for the following conjectures 

X2 '° = 1 ' X2 < 2 = (k + l)(k + 2) (25) 
and moderately strong support for the additional conjecture 

16 , , 

X4 ' 2 = (* + !)(* + 2) • (26) 

The denominator (k+l)(k+2) of these expressions is in accord with heuristic arguments; operators 
arising in the unitary coset conformal field theories (|12|) describing these critical points indeed 
have dimensions with this denominator. 



4 The critical phase 

The phase diagram for A < A c is different for the height and loop models. To understand this, let 
us first study the region around A = — 1, perturbing away from the decoupled antiferromagnetic 
self-dual critical point. The situation is much subtler than at A = 1, because the models for 
A < —1/d are not unitary. We argue here that the perturbation around A = —1 is relevant for 
the height model, but irrelevant for the loop model. 

Recall that the operator eiE2 responsible for the flow away from A = 1 has dimension (k + 
4)/(fc + 1), and so is relevant for k > 2. In the loop model, the points A = 1 and A = — 1 are 
linked by the analytic continuation k i— > — jrrp and so £162 has dimension 3k + 4 and is irrelevant 
for any k > 0. The decoupled point A = — 1 should therefore act as an RG attr actor for a whole 
range of parameter values A £ [—1, Ao), i.e., we have a critical phase with constant values of the 
critical exponents. 

The existence of a critical phase in the loop model is verified numerically in Fig. [6] for the case 
k = 3. The central charge then reads, according to (|10p . c = — =1= = —36.4, in fine agreement 
with the numerics despite of the large corrections-to-scaling effects. 

We next examine more closely the extent of the critical phase, and the nature of the transition 
at Ao. Fig. [7] shows (still for k = 3) the free energy per unit area f(L) = — i log A max for various 
strip widths L, where A max is the dominant eigenvalue of the corresponding transfer matrix. For 
L > 8 there is a jump singularity in the derivative /'(A) at some A = Xq(L), signaling a first-order 
phase transition. It is this singularity which marks the termination of the critical phase. 

The existence of this first-order phase transition can be argued analytically as well |31j . 
Indeed, at A = — 1 the k 1— ► k = — transformation implies that the dominant eigenvalue of 
the transfer matrix corresponds to a sector of non-trivial topological charge. Specifically, the 
smallest (i.e., most negative) critical exponent corresponds, for < k < 2, to a single Potts 
cluster wrapping around both the space and time directions of the lattice (and thus there are 
zero wrapping cluster boundaries), while for k > 2 it corresponds to two Potts clusters (hence 
four cluster boundaries) wrapping around the time direction. Needless to say, in the probabilistic 
regime (A > 0) the dominant eigenvalue corresponds to the identity operator, and so no cluster 
wraps around the time direction (for finite L, and in the limit M — ► 00). We conclude that for 
some Xo(L) £ (—1,0) there must be a level crossing in the transfer matrix spectrum, corresponding 
to a first-order transition. 
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Figure 6: Central charge of the loop model with k = 3, as a function of A. The pink horizontal 
line indicates the expected infinite-size behavior throughout the critical phase. 

For k = 3 we have numerically determined 

A (L) = -0.72793, -0.62225, -0.56084, -0.51985, -0.49002 for L = 8, 10, 12, 14, 16. 

In the L — > oo limit this extrapolates to A = —0.25 ± 0.05, in conspicuous agreement with 
A c ~ —0.2212 from (jlip . Similar agreements are found for other values of k. 

We are therefore led to conjecture that Ao = A c for any k. This means that at the critical 
point A c the loop model stands at a continuous first-order phase transition [32], while the height 
model stands at a conventional second-order transition. The two transitions are in the same 
universality class. This scenario is closely analogous to the one found for the antiferromagnetic 
transition in a single Potts model |33|. [34] . 

Turning now to the height model, a perturbation away from A = —1 in the Zj. parafermion 
theory ([9]). As with the perturbation around A = — 1, the perturbing operator here must be odd 
under the individual dualities of the models, and even under the combination. In the unitary 
Hamiltonian formulation, this corresponds to the product of energy operators in the parafermion 
theories, which is relevant with dimension 4/(& + 2). Accordingly, we expect a different scenario 
in which the critical phase does not exist. Numerical transfer matrix diagonalization reveals that 
the free energies f(L) are smooth functions all the way down to A = —1. For A > Xo(L) they are 
identically equal to those of the loop model. Exactly at A = — 1 there is a jump singularity in 
the derivative /'(A), as is consistent with the overall A h- > 1/A symmetry of the model. [At A = 1 
the Perron- Frobenius theorem prevents such a singularity from occurring.] In between A = — 1 
and A = A c the behavior is non-critical, the effective central charge being zero. 
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Figure 7: Free energies f(L) of the loop model with k = 3, as functions of A. Level crossings are 
visible for strip widths L > 8. 

5 Dilute loops 

The loop models we have studied thus far are completely packed: every site of the lattice is 
covered by strands. It is natural to relax this constraint, and consider "dilute" loops. We show 
in this section that the critical point and phase still persist for fairly large amounts of dilution. 




Figure 8: The five additional configurations at each vertex in the dilute loop model. 

The dilute model we discuss here requires that loops of the two species be on the same links: 
they are "glued" together. This does not reduce to a single loop model, because when two loops 
meet at the same vertex, the four possibilities in Fig. [3] are still possible. For reasons we will 
discuss below, the physics does not seem to be changed by whether or not one allows vertices 
with the pair of loops going straight through. Thus for simplicity, we still require that loops 
turn by 90 degrees at each vertex, so that there are nine allowed configurations at each vertex. 
The five new possibilities are displayed in Fig. [HJ in addition to the empty vertex, there are four 
vertices where two links touching each vertex are empty. 
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Let the weight of the vertices in Fig. [8] (from left to right) be 



u>i = 1 and UJ2 = uj-t, = W4 = W5 = K 



(27) 



and those of the vertices in Fig. [3] (from left to right) be 



uj§ = u)7 = K 2 (3 and us = u>9 = K 2 (3\ 



(28) 



Then K roles as the fugacity per monomer, (3 is a contact interaction, and A occurs when different 
loop species cross as before. We expect the contact interaction to be irrelevant, and so we set 
(5 = 1 in what follows. The dilute coupled loop models therefore have partition function 



where iV m is the total number of monomers, and the sum is over dilute loops with the vertices 
described above. The completely packed model is recovered in the limit K — > 00. 

A critical point in this dilute loop model was discussed at length in [T7] . There several indirect 
arguments were given as to why this dilute loop model will have a critical point described precisely 
by the same coset conformal field theory (|12[) as we showed above described the critical point at 
A = A c in the completely packed limit. In this section we review these arguments, use numerics 
to provide further strong support for them, and map out the full phase diagram of the dilute 
loop model. 

The arguments presented in [17] were the reverse of those presented above. In [17J, critical 
theories known to describe the coset models (|12p were shown to be closely related to this dilute 
loop model. One argument involves an integrable lattice height model |35[ [36] which has a critical 
point whose continuum limit is described by the coset (|12p . It was shown that in a certain off- 
critical limit, the height model reduces precisely to the dilute loop model. An integrable line of 
couplings connects this limit to the critical point, so that the corner-transfer matrix eigenvalues 
and hence the local height probabilities can be computed along this line. The configurations 
which dominate this computation are precisely those of the dilute loop model. Moreover, the 
remaining configurations can be described by a more general dilute loop model, where multiple 
strands of each species are allowed on a link. The fact that the configurations with multiply 
occupied links do not dominate the corner-transfer matrix computation is a strong sign that at 
the critical point, allowing or disallowing such configurations is irrelevant. Thus the dilute loop 
model discussed above should still have a critical point described by (fT2|) . 

Another argument briefly mentioned in [17] comes from studying the S matrix of the integrable 
field theory at and near this critical point. A compelling picture of the relation between loop 
models and integrable S matrices was presented in [37 1 . There it was argued that one could think 
of the worldlines of particles in the integrable 1 + Id S matrix description as dilute loops in the 
corresponding classical 2d lattice model. This observation was made precise in [38J, where it was 
shown that the S matrices for the perturbation of the minimal models can be expressed in 
terms of the Temperley-Lieb generators, which in turn have the loop representation described 
in section ([2]). For the S matrix S p for the minimal model of central charge c = 1 — 6/[(p + 
l)(p + 2)], the weight of closed loops in this representation is 2cos(-7r/(p + 1)) (note the shift in 
the denominator from the weight of the loops in the completely packed lattice model) . It is very 
natural to think of these world lines as domain walls in the 2d classical picture, and this S matrix 
indeed describes the loop model in its dilute phase. Off the critical point, these domain walls 
have an effective weight per unit length. 




(29) 



C M 
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The coupled loop models can be discussed in the same 5 matrix picture. In fact, it is no more 
difficult to treat the more general case 

SU(2) k x SU{2) 1 

SU(2) k+l m 

which will be discussed in more detail in the next section. This conformal field theory is integrable 
when perturbed by the (l,l;adjoint) operator (the analog of ^3 in the minimal models) of 
dimension 2(1 — l/(k + I + 2)). This integrable field theory describes the lattice height model 
of [35] near the critical point, along the line of couplings mentioned above. The 5 matrices for 
this field theory have the interesting property that they decompose into the tensor product of 5 
matrices of the minimal models S k X Si |39j . Each of S k and Si involves a single Temperley-Lieb 
algebra, so their tensor product can be written in terms of tensor products of Temperley-Lieb 
generators, namely 1(8)1, e(g) /, 1®/, and e(g) 1. The loop representation of these is exactly that 
pictured in ([3]) above! Thus, the heuristic interpretation of this S matrix is that it describes a 
loop model exactly like ours. Since the 5 matrix is a tensor product, each particle has structure 
in both Temperley-Lieb algebras, so its worldlines must be described as doubled lines, and the 
lines of the different copies can cross at a vertex. 

Still more evidence that the dilute loop model has a critical point comes from another solvable 
model [lUJ closely related to our dilute loop model. The models in [3D] are based on the "dilute 
BMW" algebra, generalizing those of [24J with weights (I20p to allow heights on adjacent sites 
to be the same. In the graphical representation, this amounts to relaxing the complete-packing 
requirement. As we have described above in section the 50(4) BMW algebra decomposes into 
the direct product of two Temperley-Lieb algebras. Thus the dilute 50(4) BMW algebra in its 
loop representations describes the same dilute loop model. To make it integrable, one must tune 
various couplings, including that of the vertex where the loops go straight across. Unfortunately, 
the corner transfer matrix computation of these models does not seem to have been done, so the 
CFT describing this critical point has not yet been identified. However, by using the inversion 
relation method, one can easily compute a critical exponent at each of the two critical points in 
this model. At the critical point where the weights favor dilute loops, the critical exponent found 
corresponds to an operator of dimension 3/(k + 2), which indeed is present in the coset model 
(I12p . Unfortunately, at the isotropic point some Boltzmann weights are negative (although the 
corresponding 1 + 1 dimensional quantum Hamiltonian is hermitian), so it is not possible to make 
a precise comparison with our dilute loop model. 

All this is compelling evidence that the dilute loop has a critical point with the same conformal 
field theory description as that of the completely packed model at A = A c . Minor modifications 
in the definition of the the loop model, such as allowing multiple strands on a link, or loops to 
go straight through a vertex, do not change the conclusion. 

There are no marginal operators in this coset conformal field theory describing the dilute and 
completely packed critical point. Thus there should be a flow in the loop model connecting the 
two points. 

Several compelling pieces of evidence imply that the completely packed critical point is stable, 
i.e., including dilution is an irrelevant perturbation. It has long been known that dilution is 
irrelevant in the context of a single Potts model. An initial Coulomb gas-type (and thus non- 
rigorous) argument [41] was subsequently placed on firm ground through the integrability analysis 
of the most general model of self-avoiding loops defined on the square lattice that is compatible 
with an 0(n) symmetry [421 143] . For a given value of n £ (—2, 2), this model admits five distinct 
integrable points with isotropic weights at the vertices. One is the completely packed model 
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Figure 9: Central charge c as a function of K in the dilute model. The parameter k = 3, 
corresponding to c = |jj = 1.35 (shown as a pink horizontal line) at A c = — \/2sin(^) ~ —0.2212. 
Curves of the same color correspond to several different values of A = —0.1, —0.2, —0.3, —0.4, —0.5. 
As A decreases the curves are shifted towards the right. 

which is equivalent to the Q = n 2 state Potts model. Two other solutions give the so-called dense 
and dilute solutions for the O(n) model, and fall in the same universality class |44[ 145] as those 
found previously [46] for the honeycomb lattice [47]. Crucially, the dense solution is in the same 
universality class as the Potts model solution, and so we learn that a deviation from complete 
packing is an irrelevant perturbation. [The remaining two solutions are multicritical extensions 
of the dense and dilute solutions.] 

Returning to our case of two coupled models, we thus have that at the decoupled critical 
points at A = ±1, the flow in K is toward the completely packed line (K = oo). It is natural to 
assume that this persists for all A, so that in particular dilution is irrelevant both throughout the 
critical phase discussed in section [4] and near the critical point discussed in section [3l 

To complete this picture we turn to numerical transfer matrix calculations. We have been able 
to study the dilute model (f29l) for periodic strips of width up to L = 140 Results for the central 
charge c as a function of K, but fixed A, are shown in Fig. [9j still for the case k = 3. Similar 
results hold true for other parameter values k > 2. The various curves correspond to different 
choices of A. The reader may train his/her eye to view this as a three-dimensional plot of c as a 
function of the two parameters (K, A). Applying Zamolodchikov's c-theorem — which states that 

2 Note that the exclusion of straight-going vertices from Fig.[8]is very convenient, since it reduces the dimension 
of the Hilbert space underlying the transfer matrix construction. Indeed, each element of the space is a pair (one 
for each loop species) of non-perfect matchings of L points, subject to the usual constraint of planarity with respect 
to each loop species. In the absence of straight-going vertices, each element in the matching is a pair of points on 
sublattices of opposite parity. With straight-going vertices allowed, there is no such parity constraint. 
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Figure 10: Schematic phase diagram for the dilute model with k > 2. 

RG flows go downhill — reveals that for a large enough K the flow is towards (K, A) = (oo, A c ), i.e., 
the critical point of the completely packed model. By "large enough" is meant that K > K C (X), 
where K c (\) is the locus of the global maximum in c. On the other side of the phase boundary, 
K < K C (X), the flow is towards a non-critical phase, e.g. the trivial fixed point where all loops 
are absent, or the (K,X) = (oo,0) single Potts model point S. It is readily verified that dilution 
is also irrelevant throughout the critical phase. Simulations on the line A = show that for k > 2 
dilution decreases the effective c, so the flow is towards K = 0. 

We therefore arrive at the schematic phase diagram depicted in Fig. [TUJ The RG flows are 
shown as arrows. A continuous first-order phase transition occurs along the dashed line. The 
points S (at A = 0), C (at A = A c ) and Da (at A = — 1) are the same as in the completely 
packed loop model (see Fig. Hj). The point M is a putative multicritical point, unstable in both 
parameters. For this precise lattice model, we do not know if there is an dilute integrable point 
like in the lattice models in the same universality class discussed above. If there is, it would be 
somewhere along the dashed line. Note that there should be an RG separatrix joining M and S, 
as shown. 

We have not been able to identify M exactly using our numerical scheme. Because of the 
A \—* 1/A duality, it is natural to assume that M is located on the decoupling line A = — 1, 
but the interference with the first-order transition line impedes the numerical verification of this 
expectation. Also, determining where the various phase boundaries end, or coalesce, is not an 
easy matter. We have however verified that the dilute critical phase extending from C is also 
stable towards a perturbation in the hitherto ignored parameter (3 of (|28p . 

If we generalize the coupled dilute model to include all the interactions present in the most 
general O(n) symmetric loop model referred to above |42| I43j we certainly have three candidate 
multicritical points along the decoupling line A = — 1 (viz. dilute, multicritical dense, and multi- 
critical dilute). Our constraints on the parameter space are however likely to make at least two 
of these points disappear. On the other hand, moving away from A = — 1 could in principle lead 
to a different kind of multicriticality. 
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6 The self-dual critical point for Q\ ^ Q2 



At the beginning of section (H we reviewed the arguments of [17] indicating that the dilute coupled 
loop models had a critical point. Both the S matrix and the corner transfer matrix arguments 
apply to a more general class of coset models, those built on the cosets (|30p . The central charge 
of this critical point is 

m{A + k + i) 

C (2 + fc)(2 + l)(2 + fc + Q 1 ' 

In the loop language the parameters k and I parameterize the weights for the two different species 
of loops: 

/ — f ^ \ 1 — ( 7r 

Y Qk = 2 cos I I and \J Qi = 2 cos 



k + 2) v ^' U + 2 

These correspond simply to replacing (|29|) by 

z = Y,Yl ^ Nm QT t2 QT l2 ^ nx (32) 

c M 

so that we are now coupling a Q^-state and a Q/-state Potts model. 

It is thus natural to assume that this loop model with k 7^ I will have the same phase diagram 
as that established above for k = I. In particular, there should be a non-trivial critical point in 
the completely packed model at some negative value of A, described in the continuum limit by 
the coset conformal field theory (|30p . 

We have checked this numerically for several integer values of the parameters k and Based 
on the maxima of the central charge (cf. Fig. [5]) , we find in particular 

c = 1.4675(4) at A c = -0.2908(7) for (k, I) = (3, 4) 

c = 1.5431(3) at A c = -0.3376(2) for (k, I) = (3, 5) (33) 

This agrees nicely with (i31~|) which predicts c = || = 1.46667 for (k,l) = (3,4) and c = || = 
1.54286 for (k, I) = (3, 5). We have unfortunately not been able to provide a convincing conjecture 
for A c to generalize the exact result (|lip . 
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